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ABSTRACT 

The production rate of compact objects, i.e. neutron stars (NS) and black holes (BH), in active 
galactic nuclei (AGN) and quasars (QSO), where the frequent supernova explosion is used to explain 
the high metallicity, is very high due to the interaction between the accretion disk and main sequence 
stars in the nucleus of the quasar. The compact object-red giant star (RG) binaries can be easily formed 
due to the large captured cross-section of the red giant stars. The (NS/BH, NS/BH) binary can be 
formed after the supernova explosion of the (NS/BH, RG) binary. Intense transient gamma-ray emission 
(gamma-ray burst) and gravitational radiation can result from the merger of these two compact objects. 
Collision between helium core (He) of RG and black hole may also take place and may also result in long 
duration gamma-ray bursts but no gravitational waves. We estimate that the merger rate of (NS/BH, 
NS/BH) binaries and (He, BH) is proportional to the metal abundance (7^7) and can be as high as 
1CT 3 /0.01) per year per AGN/QSO. 

Subject headings: quasar: emission line - star: neutron star/black hole - gamma-ray: burst - 
gravitational radiation 



1. INTRODUCTION 

The basic physical scenario of active galactic nuclei and 
quasars is generally believed to be a supermassive black 
hole surrounded by accretion disk which releases the ob- 
served huge body of energy (Rees 1984). However in this 
model one prominent problem is the fuel. Syer, Clark, 
& Rees (1991) investigate the possibility of accretion disk 
capturing stars in dense star cluster of nucleus (within 
lpc) in order to remove the accretion fuel problem. This 
interesting interaction has been employed to explain the 
high metallicity phenomena in active galactic nuclei and 
quasars by Artymowicz, Lin & Wampler (1993), Zurek 
et al (1994), especially for broad absorption line (BAL) 
quasars with extremely high metallicity (sometimes sev- 
eral hundred times that of solar abundance, Shields 1996). 
It is interesting to note that this process also result in 
the formation and coalescence of compact objects in the 
central energy house of active galactic nuclei and quasars 
which may be sources of frequent 7-ray and gravitational 
wave bursts. 

Gravitational wave could soon be detected directly even 
emitted from a source located at a cosmological distance 
because several ground-based gravitational wave experi- 
ments including the Laser Interferometer Gravitational- 
Wave Observatory (LIGO) (Abramovici et al. 1992), 
VIRGO, TAMA300 and GEO600 are under construction. 
The best understood of all gravitational-wave sources are 
coalescing, compact binaries composed of neutron stars 
and black holes. These NS/NS, NS/BH and BH/BH bi- 
naries may give information about the equation of state 
(EOS) of nuclear matter at high densities (Shibata, Naka- 



mura & Oohara 1992, 1993; Rasio & Shapiro 1994; Zhuge, 
Centrella & McMillan 1994; Davies et al. 1994; Ruffert, 
Janka & Schafer 1996; Cheng & Dai 1998). 

In addition the observations of afterglow of Gamma- 
ray Burst (GRB) lend strong supports to the cosmological 
origin of GRB (Paczyhski 1986, Mao & Paczyhski 1992, 
Meegan et al 1992, Piran 1992) (see a concise review by 
Greiner 1998 for the afterglow, and Piran 1998). Despite 
the effort of searching more than 20 years to look for the 
counterparts of GRB at other wave bands, no known pop- 
ulation of objects in the Universe have been identified with 
any GRB sources (e.g. Fenimore et al. 1993, Fishman & 
Meegan 1995, Band & Hartmann 1998). Two main rea- 
sons to cause this situation. First, the position accuracy of 
GRB detected by BASTE is so low that a large number of 
known objects are located in its error box. Secondly, the 
burst duration is so short that it is very difficult to use 
other means to reduce the error box. Although the Bep- 
poSAX has observed a dozen of GRBs' afterglows, one of 
these sources can reduce its error box down to ~ one arc 
second. But still no confirmed counterparts are identified, 
even for the well-observed, GBR970228 and GRB970508. 
On the other hand, the observed isotropic and inhomoge- 
neous distribution of GRBs detected by the BASTE in- 
strument on board the Compton Gamma-Ray Observa- 
tory (CGRO) is best explained by a source population at 
cosmological distances. 

Most cosmological GRB models involve compact ob- 
jects. For examples, (1) neutron star - neutron star/black 
hole merger (e.g. Paczyhski 1986; Eichler et al 1989; Der- 
mer 1992; Narayan, Piran & Shemi 1991; Mao & Paczyhski 
1992); (2) collision between black hole and Helium core or 
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white dwarf by Fryer & Woosley (1998), Fryer et al (1998); 
(3)the phase transition from neutron star to strange star 
(Cheng & Dai 1996); (4) the transition of normal nu- 
clear matter to matter with pion condensation in neu- 
tron star (e.g. Hacnscl et al 1990, Muto et al 1993); (5) 
rapidly rotating neutron star with extremely large mag- 
netic field (Usov 1992, Duncan & Thompson 1992); (6) a 
failed supernova type lb (Woosley 1993), or microquasar 
(Paczyhski 1998) etc. Therefore, if we want to find the 
hosts of GRB sources, it is logical to look for cosmological 
sources with a lot of compact objects. From the stellar 
evolution point of view, the progenitors of compact ob- 
jects, e.g. neutron stars, are massive stars (approximately 
5 ~ 20M Q ). ^From the observations of afterglows of GRBs 
it has been deduced by Paczyhski (1998) that GRB may 
associate with the star formation region. We note that 
most of the masses of the massive stars, consisting of a lot 
of heavy elements, will be ejected during the supernova ex- 
plosion and a metal-rich region will be formed. Naturally, 
galaxies with high abundance of compact object should be 
metal-rich. Therefore it is logical to look for cosmologi- 
cal sources with high metal abundance. In order to ob- 
serve the sources located at cosmological distance, these 
sources must be luminous. It should be noted that most 
AGNs/QSOs are metal-rich (Artymowicz, Lin & Wamplcr 
1993, Hamann et al 1997). This may be the real reason 
why GRBs may be related to AGNs and QSOs as reported 
by Schartel, Andcrnach & Grciner (1996), and Burenin et 
al (1998) who claim there are possible association of 7- 
ray burst with radio quiet quasar, and with active galactic 
nuclei, respectively. However the physical scenario of this 
association remains open if it is true, at least, some of GRB 
are related with AGNs/QSOs. It has been suggested that 
GRB may be produced in a star formation region which 
leads to expectation of the burst rate of GRB proportional 
to the formation rate of star in the cosmological galaxies 
(Paczyhski 1998). Because the lifetime of AGN/QSO is 
rather short (generally less than 10 8 yr), the possible host 
objects of GRB are mostly ones related with galaxies. Al- 
ternatively here we suggest that the interaction between 
accretion disk and stars in dense cluster may lead to the 
formation of a large number compact objects in the central 
engine of active galactic nuclei. The formation and inter- 
actions among the compact objects in the nuclei might 
associate with GRB. Therefore active galactic nuclei are 
one kind of its host objects if GRBs are indeed associated 
with compact objects. 

2. THE MODEL 

The mechanism for the enhancement of metal-rich ele- 
ments has been suggested by a way of evolution of main 
sequence stars, which are captured by the accretion disk 
surrounding the supermassive black hole in the center of 
quasar (for a review cf. Lin & Papaloizou 1996). The 
captured stars would accrete material from the disk and 
increase their masses up to 50 solar masses in less than 
10 4 years. These massive stars evolve rapidly toward the 
supernova stage and the ejecta of the supernovae provide 
enough heavy element to produce the features suggested 
by observations. The metal abundance expressed in terms 
of the ratio between NV and CIV is related to the to- 
tal number of captured main sequence stars Nms during 
the active phase of quasar as follows (Artymowicz, Lin & 



Wampler 1993; Zurek et al 1994) 
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However this mechanism results in other implications, such 
as formation of neutron stars/black holes in the vicinity of 
the accretion disk. As we will argue in the subsequent 
section, it exists an efficient way to form compact object 
- compact object (i.e. NS/BH, NS/BH) binaries via the 
formation of compact object - red giant star (NS/BH, RG) 
binaries. The (NS/BH, NS/BH) binary will be formed af- 
ter the explosion of the red giant star. In such a scheme 
the merger rate of (NS/BH, NS/BH) binary approximately 
equals to the capture rate of main sequence star, therefore 
it is proportional to the metal abundance. If the merger 
of (NS/BH, NS/BH) binary or Helium core of RG and 
black hole is one of the possible mechanisms of GRBs, 
then active galactic nuclei and quasars should associate 
with GRBs. In fact the coalescence of two compact ob- 
jects also produces intense transient gravitational waves 
which should be detected by LIGO/ VIRGO in lGpc ra- 
dius. 

2.1. Evolution of the captured stars in the accretion disk 

We assume that the nucleus of the quasar consists of a 
massive black hole surrounded by an accretion disk, which 
are enclosed by a star cluster. The structure of disk is 
described by the standard model in which the local height 
and the surface density of the accretion disk are given by 
(e.g. Frank, et al 1992) 

4") = 1.6 x lO-^- 1 / 10 ™ 3 / 2 ^- 1710 ^ 8 , (2) 
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and 



E = 1.7 x 10 7 a- 4 / 5 m 7 /io Mg i/5 r -3/4 (g cm _ 2); (g) 



respectively, and the inward drift velocity due to the out- 
ward transportation of the angular momentum in disk 
reads 

Vl = 4.1 x 10 4 a 4/5 ™ 3 / 10 M 8 - 1/5 r7 1/4 (cm s" 1 ), (4) 
or the time scale of inward drift is 

r r = 4.0 x Wa-^rn-^Mf 5 [^f" (yrs), (5) 

where m is dimcnsionless accretion rate scaled with the 
Eddington limit, a is the viscosity coefficient, r<j = 
is the radius of disk normalized by Schwartzchild radius 
R s = GA ^F ii , and M 8 denotes the mass of central black 
hole (M BH ) in units of 1O 8 M . 

The main sequence stars in the star cluster will be cap- 
tured by the disk and the number of the captured main 
sequence stars in the annular R— R + dR is estimated as 
(Artymowicz, Lin & Wampler 1993), 



dN = m*G 2 (ERdR) 



327T 1 / 2 O AT 



0.45 



^0 \ / AT \ / <jq 
K M Q J Vl0 7 pc~ 3 / V10 8 yr/ UOOkms" 1 

a -4/ 5rh 7/10 M ll/5 r l/4 rfrd 



(6) 



3 



where AT is the active period of quasar, vq is the number 
density of star cluster at 1 pc in the range of about 10 7 /pc 3 
(Blandford 1991), oo is the dispersion velocity of star in 
the cluster, and Cd ~ 6. Integrating over the disk, there 
are ~ 10 5 main sequence stars captured by disk in the ac- 
tive phase of quasar or the captured rate is ~ 10 -3 yr~ 1 . 
The captured stars will coorotate with the medium in the 
disk (Syer et al 1991), and accrete gas from the disk with 
Bondi rate until a gap appears in disk with condition that 
the star's Roche radius exceeds the disk scale-height h (Lin 
& Papaloizou, 1996). The timescale of Bondi accretion is 

r a = ITOxIOV/^-Viom^ (^ 3/ 4 ( ™ V 1 



10 5 / \M, 



(yrs), 

(7) 

and the maxium accreted mass of the captured star is lim- 
ited by 

^-^a^V^Mr^f 8 - (8) 

This condition is also coincided with that the Bondi ra- 
dius (r-Q = Gm^/c 2 ) does not exceed the local height h. 
With the typical parameters of the quasars and a = 1, 
the maximum accreted mass of the captured stars is in 
the range of 10 ~ 20 M Q (but the maximum mass of the 
captured stars could be much higher than 20 M© if a is 
much less than unity). The star in this mass range will 
evolve off their main sequence quickly and likely become 
neutron stars plus small fraction of black holes (see e.g. 
Shapiro & Teukolsky 1983). For more massive stars (m > 
20 Mq) they could become neutron stars or black holes 
(Timmes, Woosley & Weaver 1996). The simple formula 
of the evolution time scale for the main sequence star is 
(Meurs & van den Heuvel 1989) 
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where the index a\ and 02 are tabulated in their table 
3. For interested cases, we have a\ — 9.3, a 2 = —2 for 
3.8 < to/Mq < 12, and a\ — 8.2, a 2 = —1 for star larger 
than 12M . Substituting parameters for lower mass case 
into Eq.9, we obtain: 



r c = 9.0 x 10 6 a 3 / 10 m- 9 / 2a M- 7/10 (^) 



-3/8 



(yrs). 
(10) 

The higher mass case will evolve faster than that given 
in Eq.(10). It has been estimated that the main sequence 
stars start to be captured by the disk at the disk radius ~1 
pc (Artymowicz, Lin & Wampler 1993) but the exact value 
of this radius is not crucial in our model. The more impor- 
tant radius R c is before which the captured stars finish the 
evolution of supernova stage. This radius must be larger 

/ \ 1/3 

than that of the tidal radius R t = ( ) r* , otherwise 

the captured stars will disrupted by the central massive 
black hole. This radius can be estimated by equating the 
radial inward drifting time scale Eq.(5) and the evolution 
time scale Eq.(ll) and is given by 

r c ^ = 4.0x 10 4 a 44 / 65 m- 6 / 65 M 8 - 76/65 . (11) 

R s 



It appears that this radius is always larger than that of Rt 
for typical quasar parameters. 

2.2. Formation of binaries 

In the early stage, almost all of the captured stars can 
evolve to neutron stars/black holes with mass m ns which 
will be ejected from the disk after supernova explosion. 
But these neutron stars/black holes can only shoot up to 
a scale height 

v 2 R 2 

- < 12 > 

where Voc ~ 300 km/s is the neutron star kick velocity 
produced by supernova explosion. Taking the typical pa- 
rameters, h ns ~ 10 16 cm. The compact stars will be os- 
cillating up and down crossing the disk. Artymowicz, Lin 
& Wampler (1993) also note the possibility of retrapping 
of compact stars by disk, but the following case is more 
interesting. Before evolving to compact star, the captured 
star will go through the red giant phase with radius (-Rrg) 
and mass Mrg- During the close encounter with a neutron 
star, the red giant star would undergo substantial tidal de- 
formation at the cost of a part of the relative kinetic energy 
of the orbit. Such a tidal process can eventually dissipate 
the total positive energy of the initial unbound orbit via 
oscillations and heating, and a (NS/BH, RG) binary sys- 
tem will be created. One might argue that the velocity 
dispersion of NS/BH is much larger than the escape ve- 
locity of RG at its surface (about 70km/s). However, it 
is interesting to note the following situation. The surface 
density of RG Erg ~ Mrg/.R rg ~ 10 7 g/cm 2 is much 
larger than that of disk ~ 10 3 g/cm 2 at O.lpc. This leads 
to an efficient dissipation of NS /BH kinetic energy in the 
process of encountering RG's envelope. The probability 
of NS/BH encountering the dense envelope of RG can be 
easily estimated 
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where Sr G is the total surface area of all RG, and Sdisk is 
the area of disk capturing main sequence star. The drag 
force acting on NS/BH by the dense envelope of RG reads 
(Artymowicz et al 1993) 



F d 



(14) 



After N times of encountering the disk, the kinetic energy 
of NS/BH will be reduced significantly so that NS/BH can 
be captured by RG through tidal process, 
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(15) 



and the capture time scale of NS /BH by the RG envelope 
is about 



T„ S .RG = ^V— ~ 10 5 yr, 

Voo 



(16) 



4 



which is slightly shorter than the life time of RG, there- 
fore a large fraction of RG can be captured by the compact 
objects in the vicinity of the disk to form (NS/BH, RG) 
binaries if we take into account the interaction between 
NS/BH and the envelope of RG. Let us consider the tidal 
capture rate in more detail. The distance D for the tidal 
capture to form a (NS/BH, RG) binary can be estimated 
as (Bhattacharya & van den Heuvel 1991) 



D « 1.75i? RG 



+ M RG R q 



Mi 



RG 



Mr, 



R 



RG 



/ 50Km/s 



where v ns is the reduced velocity of the compact star. The 
cross section of this capture process is given by 



irD 2 



1 + 



2G(m ns + M RG ) 



(18) 



and is ~ 3 x 10 28 cm 2 for typical parameters. The for- 
mation rate of (NS/BH, RG) binary thus can be written 
as 



dN 



NR. 



dt 



Nn 

V m 



= i — ) [j^ ) ::/' V, V i;i ;. (1!)) 



where 7V ns and V ns ~ h ns R 2 (N RG and Vrg ) are the num- 
ber and occupied volume of compact stars (red giant stars) 
respectively, and constant k — u n sc/^ns- It should be 
noted that N ns and Arg are time-dependent because the 
presence of (NS/BH, RG) binary formation results from 
the tidal interaction between compact stars and red giant 
stars. The time-dependent numbers of compact stars and 
red giant stars obey the following equations 



dN RG 



dt 



r - kN DS N RG , 



dN a; 
dt 



N RG 
trg 



kN ns N RG , 



(20) 



(21) 



where T is the capture rate of main sequence stars by disk, 
and trg represents the evolution time scale from red giant 
star to NS/BH (about a few times of 10 5 years). These 
non-linear equations can be analytically resolved in some 
interesting stages. 

(1) When t < jt, the production of red giant stars dom- 
inates, we have 

A RG » Ft, and A ns w ( — ) t 2 . (22) 

2) In steady state, the numbers of compact stars and 
red giant stars reach constants, i.e. 



Arg(oo) = trcR and A ns (oo) 



1 



fen 



RG 



(23) 



It is clear that the (NS, RG) binary formation rate is a 
constant, i.e. = T. 

2.3. Coalescence of (NS/BH, NS/BH) and (BH, He) 

binaries 

The evolution of massive binaries leads to two possibili- 
ties: (1) Merging of Helium core of RG and black hole, (2) 



(NS/BH, NS/BH) binary formations and merging. The 
first case occurs for very close binaries in which the black 
hole companion encounters with the Helium core before 
the secondary evolves into compact object. This process 
is recently suggested as a GRB mechanism by Fryer & 
Woosley (1998). The second corresponds to the case which 
the secondary evolves into compact star faster than that 
of spiral-in process. Then the (NS/BH, NS/BH) binaries 
are formed and eventually merge to emit intense transient 
gamma-rays and graviational waves. The estimations of 
above processes are given below. 

The captured compact star may spiral-in the red giant 
star with the core mass M coro and envelope mass M en . 
Many detail calculations have been done (e.g. Boden- 
heimer & Taam 1984; Taam & Bodenheimer 1989, Taamr 
et al 1997). A rough estimation of final separation of the 
binary after spiral-in can be made by comparing the bind- 
ing energy of envelope with the difference in total energy of 
the binary before and after spiral-in (Webbink et al 1983, 
Verbunt 1993) 
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Di 



M r , 



M RG 



1 + 



2Dj M e: 
a*Ai?L m n 



(24) 



where Di, and Df arc the initial and final separation 
respectively, Rl is the Roche radius, which is given by 
(Eggleton 1983) 



Rl 
D, 



0.49 



0.6 + g- 2 / 3 ln(l + g 1 /3)' 



(25) 



where q is the mass ratio of (NS/BH, RG) binary sys- 
tem, a* is the efficiency with wich released binary energy 
is used to lift the envelope; A a factor depending on the 
mass distribution in the envelope of RG. Their product is 
often taken to be 0.20. The core mass of red giant star for 
M RG > 7M Q is given by (Iben 1993), 



M 



M \M Q J 



1.57 



(26) 



And the radius of core is given by (Paczyhski & Kozlowski 
1972; Pols & Marinus 1994; Wijers 1998) 



-RrT-^vwr 



0.6 



(27) 



If Df is smaller than the sum R core + R ns w R CO ic, after 
spiral-in the compact star merges with the core of the red 
giant star. We have calculated some typical cases listed 
in Table 1. It can be found that the RG massive than 
18Mq will merger with the Helium core before they form 
a (NS/BH, NS/BH) binary. Even the companion star is a 
neutron star, it may accrete efficiently material from the 
red giant to become a black hole during the spiral-in pro- 
cess (Chevalier 1996; Fryer & Woosley 1998). Therefore the 
merger of (NS/BH, He) binary can become 7-ray burst as 
described by Fryer & Woosley (1998). 

There are two possible ways to destroy the formed 
(NS/RG, RG) binaries: 1) the combining gravitation of 
central black hole and the dense star cluster, 2) the third 
encounter from the cluster. It has been shown by Saslaw 
(1985) that the chance of the later way is rather small. 
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We can now estimate the action of the first way. If the 
tidal force from the gravitation of central massive black 
hole and dense star cluster exceeds the interacting force 
between the two components, then the (NS/BH, RG) will 
be disrupted by the tidal force. The tidal force reads 



-Ftid — 



G(M BH + M sc )(m c + M RG ) Di 



R 2 



R 



(28) 



where M$c is the mass of dense star cluster, and m c is the 
mass of the compact object. Then the ratio between the 
tidal force and the gravitational force (Fb) between the 
two components is 



Fid = M BH + M sc Msg + m c 
Fb Mrq m c 



10" 



(29) 



which means that the formed binaries can survive in 
its surroundings. On the other hand, consequently the 
spiral- in process takes place in the binary very fast. The 
timescale of this process can be roughly estimated un- 
der the assumption that the formed (NS/BH, RG) binary 
keeps the total angular momentum constant. According 
to Verbunt (1993), this timescale reads 



1 



Mi 



RG 



' sp 



2(1 - m c /M RG ) M 



1 M RG 

2 M ' 



(30) 



where M is the rate of mass transfer onto the compact ob- 
ject. For the case of the RG as a donor in (NS/BH, RG) 
binary system, M can be estimated from Bhattacharya & 
van den Heuvel (1991), 

M = 1.4 x lO^Mj 2 ' 25 ^- 75 (M Q vr -1 ), (31) 

where M x = M RG /M , and R x = R RG /R Q . We thus 
obtain 

r sp « 35.7Mr 1M V' 75 (yr). (32) 

From the above two equations we have the typical 
timescale for the case of RG with mass 10M Q , r sp ~ 2 
yr. This simply means that the new formed (NS/BH, RG) 
binary will reach the final stable configuration Df (listed 
in the Table 1) within 2 years, and become hard from the 
soft state very rapidly. It is thus viable to form hard bi- 
nary in accretion disk. Furthermore, it has been shown 
in more detail by Magorrian & Tremaine (1999) that the 
highest tidal disruption rate 10 _4 y r_1 only take place in 
faint galaxies (L < 1O 1O L ). The tidal disruption of giant 
star could produce flare as often as every 10 5 yr. Therefore 
it is robust to believe that the present model of star-disk 
interaction works dominatcly over tidal capture by super- 
massive black hole hosting in the center of AGN/QSO in 
the active phase of AGN/QSO. 

The detailed study of orbital changes due to the ex- 
plosion of supernovae can be found in Pols & Marinus 
(1994), and Portegies Zwart & Verbunt (1996). On the 
other hand, Verbunt(1993) argued that this effect can 
be estimated by assuming that the explosion is instanta- 
neous, and the positions and velocities of the stars are the 
same after the explosion as before the explosion (Dewey & 
Cordes 1987). This implies that the distance Df between 
the compact star and the core star before the explosion 
is the periastron distance after the explosion and that the 



periastron velocity of the new orbit is the same as the or- 
bital velocity in the pre-supcrnova orbit. Therefore the 
eccentricity is given by (Verbunt 1993) 



M coro - m„ 



< I IS 



2m r 



(33) 



The mass center of binary has obtained a speed v cm be- 
cause of the mass loss and is given by 



(34) 



where v\ 



GM C , 



10 ! 



\amq) \ir & ) 



1/2 



cm s 



D f ^ y 4M0 

is the orbital velocity of the exploding component before 
explosion. This velocity will make the (NS/BH, NS/BH) 
binary shooting up to a scale height about ~ 10h ns . Fur- 
thermore, the orbit of the binary is sufficiently small so it 
should be circularized by tidal interaction and the radius 
D ns of the circular orbit is given by 



D ns = (l + e ns )D f . 



(35) 



In table I, we list some possible values of e ns and D DS by 
assuming m ns = 1.4M©. We can see that e is larger than 
unity when Mrg is larger than 15M . This means that 
after explosion the binary will be broken. For m ns = 3Mq, 
we find that e is always less than unity and Df is larger 
than R c even when Mrq equals 20M Q , in other words, 
there is no merger of massive black hole and helium core. 

There are two interesting consequences of the new born 
(NS/BH, NS/BH) binary. First, the kick velocity will 
carry the binary to leave the disk up to the distance about 
I0h ns ~ 10 17 cm. Second, since the separation of two com- 
pact stars in the (NS/BH, NS/BH) binary is ~ D ns , the 
gravitational radiation of the binary will make these two 
compact stars merge in the time scale of (Peter 1964) 



tm(D B 



12 
19 



with 



and 



Cl 



e 19 



(l-e2) 



121 
30i' 



2299 



64 G 3 M 1 M 2 (M 1 + M 2 ) 



(36) 



(37) 



(38) 



and Int(e) is the integral (see eq. 5.14 in Peter 1964). 
Therefore the compact star merger rate should be about 
r. One should not be surprised by the short time scale of 
merger since t m is proportional to D^ s while most of D ns 
are less than solar radius, and i m cx (1 — e 2 ) 7 / 2 for larger 
eccentricity. 

3. DISCUSSIONS 

In this paper we suggest that if the mctallicity of the 
metal-rich quasars result from the capture of the main se- 
quence stars by the accretion disk surrounding the central 
massive black hole of quasars and the captured stars can 
increase their masses by accretion so they can evolve off 
their main sequence branch and move toward the super- 
nova stage which provide the observed metal abundance. 
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Then the remnant stars are likely compact stars, i.e. neu- 
tron stars or black holes. We have shown that the neutron 
star/black hole density in the vicinity of the disk will be 
high enough to a population of (NS/BH, RG) bnaries whic 
can eventually evolve to become either (i) (NS/BH, NS) 
binaries if the mass of RG is less than 15M Q and merge 
in a time scale of a few million years or (ii) (NS/BH, He) 
binary. 

In a wide range of cosmological GRB models the com- 
pact star merger plays an essential role. Since the main 

sequence star capture rate T = 10~ 3 yr _1 ( N ^^ IV ) is pro- 
portional to the metal abundance and the compact ob- 
ject binary formation rate equals the main sequence star 
capture rate in our simple model. Therefore, the GBR 
burst rate -Rgrb in AGNs or QSO can be estimated as 

Rgkb ~ (39) 

where rj is the beaming factor of the 7-rays which reduces 
the observed probability of GRBs resulting from the com- 
pact star mergers (e.g. rj is 0.1-0.01 for NS/NS merger, 
Ruffert, Janka, Schafer 1996; Ruffert et al. 1997; Ruf- 
fert & Janka 1998). Such estimated GRB rate in metal- 
rich quasars is much high than that of ordinary galax- 
ies, e.g. the merger rate due to the mass exchange and 
gravitational radiation losses giving a conservative rate 
~ lO^yr" 1 (Phinney 1991, Narayan et al 1992), or at a 
higher rate 10~( 4 ~ 5 ) per year (Tutukov & Yungelson 1993, 
Portegies & Verbunt 1996) in light of the calculations of 
binary evolution. Of course, the population of AGN/QSO 
is only 1% of normal galaxies which makes the total burst 
rate for these two classes of object comparable. However 
the burst rate for individual AGN/QSO which has excep- 
tional high metal abundant could be very large. It has 
much higher possibility that the not-yet-identified host of 
GRB may an AGN/QSO with very high metallicity. 

Recently, Dokuchaev, Eroshenko & Ozernoy (1998) pro- 
pose the possibility of GRBs origin from the evolved galac- 
tic nuclei. In their model, GRBs result from the coa- 
lescence of the compact object binary which are formed 
due to the dynamical evolution of cluster and dissipation 
of gravitational radiation. Carter (1992) proposed that 
GRBs are produced by the tidal disruption by the central 
supermassive black hole. Rees (1988) even pointed out 
that the rate of tidal disruption could be 10~ 4 yr _1 which 
is significantly lower than the capture rate [see eq.(6)] and 
formation rate of (NS/BH, NS/BH) binaries, particularly 
the whole star (unless giant star) will be swallowed by 
the supermassive black hole with mass > 10 8 M Q since the 
Roche radius Rt will lie within the horizon radius of the 
black hole. Therefore it is believed in an ANG/QSO GRB 
rate due to the presently proposed process may dominate 
the tidal capture of star by the central black hole. One 



might argue that there is no AGN/QSO at the locations of 
GRB970228 and GRB970508 which are detected by Bep- 
poSAX therefore it is unlikely that GRBs are related to 
this class of object. However, in order to locate the accu- 
rate position the BeppoSAX is restricted to long duration 
GRBs. Popham, Woosley & Fryer (1998) show that only 
(BH, He) merger can produce long duration burst, other 
mergers [e.g. (NS/BH, NS/BH)] will not produce long du- 
ration burst. It is important to note that the lifetime of 
AGN is typically of 10 s years, alternatively, the present 
paper proposed that at least some of GRBs may origin 
from the active nuclei due to the interaction between the 
star in cluster and accretion disk, roughly one percent of 
the observed GRBs. In our model, only very low mass 
compact star and high mass RG binaries can form (BH, 
He) mergers. For a power law mass function, such binaries 
will be a very small fraction. It is not too surprised not to 
identify any AGN/QSO with GRB970228 and GRB970508 
which are long duration bursts. Interestingly there are two 
GRBs with small error box, in which AGN were found. It 
was reported by Drinkwater et al (1997) that the possible 
X-ray counterpart of GRB920501 is related with a Seyfert 
1 galaxy at z = 0.315. Piro et al (1998) report that the 
first X-ray location of a 7-ray burst by BeppoSAX (within 
3' radius) contains the quasar 4C 49 with z — 1.038. Fi- 
nally, we want to point out that the basic energy mecha- 
nism of gamma-ray bursts in our model are the same as 
those proposed by other authors, namely GRBs are re- 
sulting from the mergers of compact binaries, except the 
formation rate of the binaries and the location of Gamma- 
ray bursts. Since our model burst rate is proportional to 
the metal abundance (cf. eq. 39), we shall predict that 
if the GRBs can be identified with AGNs/QSOs, these 
AGNs/QSOs must have high metal abundance. 
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Table 1 The calculation of binary formations* 







Rc 


A 


Df 


e 


D ns 


im(yr) 


Type 


10.0 


2.155 


0.349 


296.021 


0.621 


0.270 


0.788 


0.804x10' 


(NS, NS/BH) 


11.0 


2.503 


0.381 


295.468 


0.613 


0.394 


0.855 


0.796xl0 7 


(NS, NS/BH) 


12.0 


2.869 


0.414 


295.003 


0.607 


0.525 


0.925 


0.635xl0 7 


(NS, NS/BH) 


13.0 


3.253 


0.446 


294.606 


0.601 


0.662 


0.999 


0.359xl0 7 


(NS, NS/BH) 


14.0 


3.655 


0.479 


294.264 


0.596 


0.805 


1.077 


0.987xl0 6 


(NS, NS/BH) 


15.0 


4.073 


0.511 


293.966 


0.592 


0.955 


1.158 


0.131xl0 5 


(NS, NS/BH) 


16.0 


4.507 


0.543 


293.704 


0.589 


> 1 






17.0 


4.957 


0.575 


293.472 


0.586 


> 1 








18.0 


5.423 


0.607 


293.265 


0.584< R c 








(He, BH) 


19.0 


5.903 


0.638 


293.079 


0.582< R c 








(He, BH) 


20.0 


6.398 


0.670 


292.911 


0.580< R c 








(He, BH) 



* All parameters are in solar units except e and t m (yr). The mass of the compact star is taken to be 1.4M Q . 



